Signatures of single site addressability in resonance fluorescence spectra 
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Pioneering methods in recent optical lattice experiments allow to focus laser beams down to a 
spot size that is comparable to the lattice constant. Inspired by this achievement, we examine the 
resonance fluorescence spectra of two-level atoms positioned in adjacent lattice sites and compare 
the case where the laser hits only one atom (single site addressing) with cases where several atoms 
are illuminated. In contrast to the case where the laser hits several atoms, the spectrum for single 
site addressing is no longer symmetric around the laser frequency. The shape of the spectrum of 
fluorescent light can therefore serve as a test for single site addressing. The effects we find can be 
attributed to a dipole-dipole interaction between the atoms due to mutual exchange of photons. 
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I. INTRODUCTION 

Over the course of the last decades, the exploration of 
the radiative properties of laser driven atomic systems 
advanced at a stunning pace in the field of quantum op- 
tics. The fluorescence light of a coherently driven two- 
level atom is a common example in most textbooks [l|, [|J . 
For a coherent laser drive, the predicted Mollow spectrum 
is a symmetrical three peak spectrum with center-side 
band separation given by the Rabi frequency and the de- 
tuning of the driving laser field Taking more than 
one atom into account, other interesting features of the 
spectrum arise due to effects of coherent and incoherent 
inter-atomic interactions [f|. The question to what ex- 
tend the fluorescence spectrum is altered in the presence 
of collective effects is thus of great interest since it con- 
tains informations about the physical setup of the atomic 
system Q. 

In this work, we clarify how the shape of the spectrum 
alters when a single atom within the atomic ensemble is 
addressed by an external driving field. Therefore we com- 
pare the spectra of the situation where a laser illuminates 
all atoms with the situation where the laser illuminates 
only one atom. The usual symmetry of the spectrum 
breaks down if the distance between neighboring atoms 
is such that their dipole-dipole interaction via mutual ex- 
change of photons is comparable to the magnitude of the 
driving strength. 

Our investigation of the resonance fluorescence spec- 
trum under the assumption that single atoms in an 
atomic ensemble can be addressed by a laser beam, is 
motivated by recent experiments with optical lattices 
[H, [9[ in a Hubbard-regime that trap neutral atoms in 
the lattice sites. High resolution imaging systems with 
an optical resolution of about 600-700 ran allow to resolve 
the fluorescent light emitted from the atoms in individ- 
ual lattice sites of a two dimensional lattice. Hence, a 
laser beam traveling on the same path just like the fluo- 
rescence light going through the imaging system, but in 



reverse direction, can be focused onto a single site or sin- 
gle atom with a full-width at half-maximum (FWHM) of 
again 600-700 nm Such single site addressing allows 
to investigate local properties of quantum many-particle 
systems 
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The origin of the effects we find for the fluorescence 
spectra of neutral atoms in optical lattices lies in the typ- 
ical inter-atomic separation. On the one hand, the atoms 
are far enough from each other to open up the possibility 
of single atom addressing by focused laser beams. On 
the other hand, the separation is small enough so that 
an inter-atomic coupling mediated by mutual exchange 
of scattered photons still influences the collective behav- 
ior. Otherwise the fluorescence spectrum would not differ 
from the well known Mollow spectrum Q . 



The remainder of the paper is organized as follows. In 
section UH we set up a master equation that describes the 
coupling of TV two-level atoms to the quantized electro- 
magnetic field. Due to small inter-atomic separations we 
have to account for an effective atom-atom coupling me- 
diated by the quantized field. The effective atom-atom 
coupling, namely the dipole-dipole interaction and the 
collective damping rate, will be discussed in detail. We 
also define our notion of single site addressing motivated 
by the experiment [Io| . As our investigations focus on 
the fluorescence spectrum of the atomic ensemble, the 
power spectrum will be introduced in section IIII1 We 
then present our numerical results in section |IV] for a 
number of up to five two-level atoms representing sin- 
gle site addressability in a ID or 2D optical lattice. In 
all cases we find a broken symmetry in the fluorescence 
spectra if only one atom is addressed by the laser field. 
We demonstrate that this effect only occurs if the dipole- 
dipole interaction is finite. In section [V] we make some 
general remarks and provide conclusions about symmet- 
ric power spectra. In section I VII we state possible ex- 
perimental applications to test single site addressability 
with resonance fluorescence measurements. Finally we 
give our conclusions in section IYIII 
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II. HAMILTONIAN AND MASTER EQUATION 

In our model we consider N identical two-level atoms 
at fixed positions r M and define the distances R Ml/ ~ 
r M — r v . We demand, however, that all atoms lie in- 
side a two-dimensional plane to which we will refer to 
as the atomic plane. The ground state of the atom \x 
is denoted by \g^) and the excited state by |e M ), where 
/i = 1,2, ...,N labels the atoms. Apart from the laser 
that illuminates them, the atoms couple to all modes of 
the quantized electromagnetic vacuum. The time evolu- 
tion of the atomic system can therefore be described by 
the master equation (l3j . 

d tP (t) = L(p(t)) = - l -[H,p(t)} 

N 

N 

+ J2 ^ fap(t)a+ - a+a~p(t) - p(t)a+a~) , (1) 

H,i/=1 

with Liouvillian L(-), where the operators cr+ = |e M )(<7 M | 
[a~ = \gfj,)(e^\) create (destroy) an excitation in the \x- 
th atom. The decay rate 7 accounts for the spontaneous 
emission from a single atom, while the collective damping 
rates Y^ v = account for the decay of the collective 
atomic states. Since the spatial separation of the atoms 
in our setup is about 1/2 of the wavelength correspond- 
ing to the atomic transition, the often used assumption 
of independent quantum environments flil Il5j for the 
individual atoms is in our study not justified. Conse- 
quently the collective decay processes associated to T^" 
are indeed important [l6|. This is also confirmed by our 
results. 

The unitary part of the time evolution of the reduced 
density matrix p(t) is governed by the Hamiltonian 

H = H a + H L + H dd , (2) 

consisting of three parts; the bare atomic part Hq, the 
laser-atom coupling Hi, and the dipole-dipole interaction 
potential H dd . In a frame that rotates with the frequency 
Wi of the driving laser the bare Hamiltonian 

h/\ N fiA N 

H ° = — E ct m = "V E (MM - ( 3 ) 

is determined by the detuning A = ujq — uji , where loq is 
the frequency of the atomic transition. The laser-atom 
coupling in turn depends on the Rabi frequency 

N 

H L = hJ2^(<J++a-). (4) 

The atoms are driven by a laser traveling perpendicular 
to the atomic plane. Hence, the laser field at each atom 
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FIG. 1. Notion of single site addressing in an optical lat- 
tice. The atoms are fixed in the lattice sites while a laser 
beam traveling perpendicular to the inter-atomic axis is fo- 
cused onto one site. The beam profile is given by a Gaussian 
envelope with a FWHM 77 as given in Eq. ([5|. The parameters 
ai and 77 match the experimental parameters of Ref. [s| . 

has the same phase. Without loss of generality we thus 
choose f2 M to be real and positive for all atoms. We con- 
sider a Gaussian beam profile focused onto the atomic 
plane and model the field amplitude at the position of 
the /j-th atom by the Rabi frequency 

2 

% = n(r fl ) = n e- Aln2 -^, (5) 

where f^o € R represents the intensity at the center of the 
beam, which defines the point of origin, with a FWHM 
of 77, see Fig. [TJ This definition is consistent with our no- 
tion of single site addressing, pictured in Fig. HQ which is 
motivated by the optical lattice experiments [8l-fl(j] with 
single atom occupation on one lattice site and where the 
central atom is addressed by a focused Gaussian laser 
beam. In our model we assume the atoms to be station- 
ary during the illumination process, i.e. we consider a 
lattice in the Mott Insulator regime. 

The coupling to the electromagnetic vacuum induces a 
dipole-dipole interaction between the atoms due to mu- 
tual exchange of photons. This dipole-dipole interaction 
arises in the derivation process of the master equation. 
Just as the collective damping rate it is a consequence 
of the fact that all atoms couple to the same quantum 
environment. The dipole-dipole interaction J^, which 
should not be confused with the interaction between two 
static dipoles, is the coherent counterpart to the incoher- 
ent collective damping rate . The operator form of 
the dipole-dipole Hamiltonian can be taken to read 

JV 

H dd = hY^ -7T- {rf a » + , ( 6 ) 

with interaction strength J^ v = J v ^ , which depends on 
only two parameters. Namely, the distance |R^jy| be- 
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tween the atoms and the angle — Z(R M „,cio) be- 
tween the inter-atomic separation vector and the dipole 
moment do of the atoms. By introducing the dimension- 



less quantity x 
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where Aq is the wavelength 



corresponding to the energy splitting cj^, wc find the fol- 
lowing form for the dipole-dipole interaction 



J^(a^,^) = |7{[cos 2 (an-l] £2 g^ 
sin(27rx^) + cos(2ttx^) 



[l-3cos 2 (aH] 



(2ttx^) 2 (2ttx» u ) 3 
and for the collective damping rate 



(7) 



r^(a^,x^) = ^7 Hi 



cos 2 (a-)] M***") 



[l- 3 cos 2 (a"")] 



2nx» u 

cos(27ra;^) sin{2itx^) 



(2ttx^) 2 (2-rrx^) 3 



(8) 



given in units of the spontaneous emission rate 7 = 

■ Both quantities decay asymptotically as ~ jjnr, 
see Fig. [3J which indicates that the spectral properties 
of many atoms separated by a large distance do not dif- 
fer from the spectral properties of a single atom. In or- 
der to find collective dynamics the inter-atomic distances 
should be at least of the order of the atomic wavelength, 
i.e. x^ v < 1. On these small length scales addressing of a 
single atom becomes challenging. We show that the col- 
lective quantities will affect the spectral response of the 
atoms and reveal whether single atom addressing takes 
place or not. 




FIG. 2. (Color online) The dipole-dipole interaction J M</ 
(solid red line) and the collective damping rate r' 1 " (blue 
dashed line) in units of 7 as a function of the inter-atomic 
separation x^ v with = f . 



the so called power spectrum. It displays the emitted flu- 
orescence intensity per energy interval and is given by the 
Fourier transform of the two-time correlation function of 
the electric field, 

S ss (r,u J )= lim ^ ^ e -^(E (_) (r,t + r)-E (+) (r,t)). 

J-oo 2n 

(11) 

By making use of Eq. (|10j) the power spectrum can be 
expressed in terms of atomic two-time correlation func- 
tions, 



III. SPECTRUM OF RESONANCE 
FLUORESCENCE 

In the following wc focus on the spectral distribution of 
the fluorescence light emitted by the atoms in the steady 
state limit. First we concentrate on the total steady state 
intensity which is given by the normally ordered one-time 
correlation function of the emitted electric field [l?} , 

I ss (r) = lim(E (_) (r,t)-E (+) (r,t)). (9) 

t— >oo 

Here, E*' + ' ) (E^ ) denotes the positive (negative) fre- 
quency part of the field operator, which is related to the 
atomic transition operator cr~ by frsj . 

e (+) M) = (10) 

47T£n C l T A 

with the retarded time t = t — - at a point r = rr in 
the far-field zone. The physical picture here is that every 
photon which is annihilated in the detection process had 
to be emitted by an atom at an earlier time t. Turning 
back to the spectral properties of the atoms we introduce 



S ss (v,lu) =S (t) lim V / dr e -''(™)- 

Uo (12) 

(a+(t + r)a-(t)) e ^<'^} , 

where } denotes the real part and 5o(r) = = 
— I ^arxirxd-o) 1 2 d eno t es a normalization factor which con- 

7r 1 47reo c-r 1 

tains the radiation properties of a dipole. The factor 
gtwir ar i ses f rom the fact that we work in a rotating 
frame with respect to the operator SuLi AH the 
atomic one- and two-time correlation functions are acces- 
sible via the master Eq. ([I]) and by usage of the quantum 
regression theorem [19( . 

To analyze the power spectrum S ss as given in Eq. (|12l) , 
it is helpful to split the expression into a coherent and an 
incoherent part. For two operators A and B the expec- 
tation value of the product AB can always be separated 
into a coherent and an incoherent part 0, [13] , 

{AB) = (A){B) + ((i - (A))(B (B))) , (13) 

respectively. It is therefore clear that the coherent part 
of the power spectrum is always proportional to a 5- 
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FIG. 3. (a) (Color online) Incoherent part of the power spectrum for 2 atoms with R21 = aie x = 0.82Aoe^, a 12 = 7r/2, J 12 ~ 
-0.097, T 12 « -O.II7, A = 7, cr = 600 nm and Qi = O.I7. The detector is positioned in the xy-plane with 6 = Z(f , R21) « 0.92. 
The solid red and dashed blue line display the spectrum under single atom addressing, i.e. f^2 ~ 0.004fii « 0, while the dotted- 
dashed green and dotted orange line consider a broad laser beam with O.I7 = fti « ^2- The dashed blue and the dotted orange 
line, however, display the case, where the collective parameters J M " and F M " have been artificially set to zero, 
(b) The degree of asymmetry as defined in Eq. (|15|) plotted as a function of the ratio for parameters as in (a) except 

for Qi = O.I7 « I J 12 \ « |T 12 | (solid red line), Qi = 0.5 7 « 5| J 12 | « 5|r 12 | (dashed blue line) and Qi = I7 « 10| J 12 | w 10|r 12 | 
(dotted black line). 



function, since 

/oo 
dre-^-^(a+(t + T))(a-(t)) 
- OO 

/OO 
dr e -i( U -^)r K _ 
-OO 

(14) 

Another feature of the power spectrum for N > 2 atoms 
is the geometry dependence arising from the exponential 
e jfc r-(r^-r„) _ p Qr a t oms w here the collective parameters 
and are negligible, this interference effect, how- 
ever, does not contribute to the interesting part of the 
spectrum, namely the incoherent part. In this case the 
expectation values of (c^"cr~) always factorizc for /1 7^ v 
and the interference terms only enter the coherent part 
of the spectrum. This results in an incoherent part of 
the spectrum which is the sum of single atom Mollow 
spectra. In cases where T^ v and J^ v are not negligible, 
the interference terms enter the incoherent part of the 
spectrum and can make it asymmetric. This effect can 
serve as a signature of single site addressing and lies at 
the center of our investigations. 

IV. BREAKDOWN OF SPECTRAL SYMMETRY 

A. Two Atoms 

To explain the basic physical mechanisms behind our 
results, it is convenient to concentrate on the spectrum of 
two atoms first. The position of atom 1 defines the point 
of origin, i.e. ri = {0,0, 0} T . Atom 2 is just positioned 



in the next lattice site at r-i = {ai, 0, 0} T . For the lattice 
constant and the FWHM of the laser beam we choose 
ai = 640nm and a = 600nm like in Ref. @. The dipole 
moments are oriented along the z-axis, do oc e z , which 
leads to a 12 = ir/2. For the D2 transition of Rubidium 
87 (transition wavelength Ao = 780nm) the collective pa- 
rameters result in T 12 sa — O.II7 and J 12 w — O.O97. 

Figure compares the case of a broad laser beam, 
i.e. fii = f^2 = O.I7, to the case of a laser beam with 
a FWHM of cr = 600nm focused onto atom 1, which 
results in fix = O.I7 w 25S12- We refer to the latter 
case as single site addressing. It is seen that the spectral 
symmetry breaks down in the case of single site address- 
ing under the effect of the dipole-dipolc interaction. For 
two atoms that are illuminated with equal intensity we 
find a spectrum of symmetric Mollow shape as shown 
in Fig. (3Ja), see dotted-dashed green line. Under single 
site addressing, however, the peak on the right hand side 
of the spectrum is suppressed which leads to an asym- 
metry, see solid red line. To illustrate that this effect 
can be attributed to the presence of the dipole-dipolc 
interaction J 12 and the collective damping rate T 12 we 
compare our results to the physically rather impossible 
situation where the collective parameters T 12 and J 12 are 
turned off artificially. As expected we find the symmetric 
Mollow spectrum; dashed blue curve for one atom being 
illuminated and dotted orange curve for two atoms. The 
latter only shows little deviations of the peak heights and 
positions from the case of 51i = H2 and J 12 7^ 0, T 12 7^ 0. 

Note that only two peaks out of the triplet (at ±A) can 
be seen because of the weak driving strength flo = O.I7, 
with Slo as defined in Eq. We choose such a weak 
Rabi frequency because the effect of symmetry breaking 
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in the spectrum is more evident if the Rabi frequency is 
of the order of the dipole-dipole interaction. In the case 
where fto J^ v the dipolc-dipolc interaction is just a 
small perturbation to the driving of a single two-level 
atom, hence the symmetric Mollow shape dominates the 
spectrum. 

Figure EJb) displays the degree of asymmetry in the 
spectra for different ratios of ^ which is equivalent to 
different FWHM of the laser beam. We consider here a 
degree of asymmetry which is defined as 

D = — ^ — (max{|S*(w) - S(-u))\ : Co = ui - lo l > 0}) , 

^max 

(15) 

normalized to the highest intensity S ma x in the spectrum, 

S m ax = ma,x{S(Co) : Co = lo - lol G (16) 

For parameters as in Fig. (3ja) the definition of the degree 
of asymmetry corresponds to the visibility of the differ- 
ence in the peak heights. One can see that the degree of 
asymmetry decreases monotonically with increasing ra- 
tios of yjj and goes exactly to zero as — > fij.. Fur- 
thermore the degree of asymmetry is highest if the Rabi 
frequency is comparable to the dipole-dipole interaction 
J 12 between the atoms. 

In conclusion, the degree of asymmetry in the spectrum 
is a signature for the degree of single site addressing, pro- 
vided that the collective parameters J^ v and Y^ v are of 
the order of the driving strength. Therefore we empha- 
size that single atom addressing in the regime of large 
inter-atomic separations is certainly not challenging but 
for small inter-atomic separations the collective atomic 
properties reveal whether single site addressing occurs in 
the system. 

Our results also give direct evidence that the dipole- 
dipole interaction J MI/ is induced by mutual exchange of 
photons between the atoms. Indeed, as already men- 
tioned the incoherent part of the spectrum only contains 
the interference of light emitted from separate atoms if 
the dipolc-dipolc interaction J^ v is finite. Yet the inter- 
ference effect is still present even if atom 2 is not illumi- 
nated by the laser at all, i.e. O2 = 0. This indicates that 
processes where atom 1 absorbs a photon from the laser 
field, emits a photon which travels to atom 2 where it is 
absorbed and emitted again must exist. Since we look at 
the spectrum in the steady state limit and the Rabi fre- 
quency Slo is comparable to J 12 these processes occur at 
a rate large enough to generate the observed interference 
in the spectral light of the atomic ensemble. Note that 
this interference is not due to coherent classical light but 
it is rather of quantum mechanical nature. 

B. Larger numbers of atoms 

If one considers the investigation of single site address- 
ing in a ID or 2D optical lattice more then one neighbor- 



ing atom should be taken into account. Since the dipole- 
dipole interaction and the collective damping rates decay 
as -jar with increasing separation x^" , one can however 
obtain estimates for larger lattices by just considering a 
limited number of lattice sites or rather atoms contribut- 
ing to the spectrum. We study the effects of single site 
addressing in a ID and 2D optical lattice up to a level 
where the atoms with the highest contribution to the 
spectrum are taken into account. These atoms are the 
ones where the dipole-dipole interaction potential J 1/J be- 
tween the /x-th atom and the addressed atom is largest 
in magnitude. 

Figure|4]displays the impact of single site addressing in 
a ID optical lattice (dashed blue curve) and a 2D optical 
lattice (dotted-dashed green curve) configuration within 
the approximation described above and compares it to 
the case of two atoms (solid red curve) and to the Mol- 
low spectrum of a single atom (black dotted curve). In 
part (a) of Fig. [4] we have chosen parameters for the 
lattice constant and the FWHM of the laser beam like 
in H, i.e. a; = 640nm and a — 600nm. For this lat- 
tice constant the magnitude of the dipolc-dipolc interac- 
tion is largest between the addressed atom positioned at 
T\ = {0,0, 0} T and the nearest neighbours in the lat- 
tice, hence the atoms positioned at r = {±a;,0,0} T 
and r = {0,±a;,0} T . The parameters A = 7 and 
Qi = O.I7 coincide with the parameters chosen in Fig- EI 
Again the detector is positioned in the xy-plane with 
= Z(r, R21) ~ 0.92. Therefore the solid red line in 
Fig. HJa) is the same as the solid red line in Fig. [3ja). 
It shows the incoherent part of the output spectrum for 
the addressed atom at ri = {0,0, 0} T in the presence of 
another atom at r = {a;,0,0} T . The dashed blue line 
displays the same situation but with a third atom placed 
at r = {-ai,0, 0} T . The dotted-dashed green line shows 
the spectrum if a fourth and a fifth atom are placed at 
r = {0,±a;,0} T . In all cases we find a broken spectral 
symmetry under single site addressing. The peak on the 
right hand side is suppressed for all atomic configura- 
tions as compared to the single atom Mollow spectrum. 
We also note that spacial symmetry of the atomic setup 
does not lead to spectral symmetry of the emitted light. 

Figure @|b) displays the same situation as Fig. QJa) 
but with lattice parameters as in the experiments of 
Ref. d, [l(|. The lattice constant in the experiment is 
equal to 532 nm. This leads to a dipole-dipole interac- 
tion of J^ v ps 0.037 between atoms positioned in two ad- 
jacent lattice sites, while atoms separated by y/2ai have 
a dipole-dipole interaction of J 1 *" ps — O.I27. Hence, in 
our numerical calculations we concentrate on the contri- 
butions from these atoms and neglect the contributions 
from the nearest neighbours of the addressed atom. Al- 
though the breaking of spectral symmetry is displayed 
in Fig. @|b) as well, we notice that the overall intensity 
of the case with a finite dipole-dipole interaction is a bit 
larger than for a single atom. Figure UJa) exhibits the 
opposite behaviour. This effect is mainly due to the dif- 
ferent positioning of the detector in Fig.fJJa) andHJb). 
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(a) (b) 




FIG. 4. (Color online) Incoherent part of the power spectrum (a), (b) with the corresponding atomic configurations in the 
lattice below, respectively (c), (d). (a) Incoherent part of the power spectrum with parameters as in Fig. [3fa) . The solid red 
curve shows the spectrum for 2 atoms under single atom addressing just as the solid red line in Fig.[3{a). The dashed blue curve 
displays the spectrum of three atoms where the third atom is placed at rs — —aie x . This represents single site addressing in a 
ID optical lattice where the contributions of the nearest neighbours of the addressed atom are taken into account. For the 2D 
lattice we place two more atoms at r4 = aie y and r$ = —aie y and find very similar features (dotted-dashed green curve). The 
Mollow spectrum of a single atom is given also for comparison (dotted black curve), (b) Incoherent part of the power spectrum 
for 2 atoms (solid red curve), 3 atoms (dashed blue curve), 5 atoms (dotted-dashed green curve) with the parameters A = 7 
and Q,o = O.57. The positions of the atoms are ri = {0, 0, 0} T , T2 = a t {\^2, 0, 0} T , r3 = ai{ — y/2, 0, 0} T , = ai{0, y/2, 0} T and 
rs = ai{0, — \/2, 0} T , as illustrated in the picture below, with a; = 532 nm and a = 700 nm as in Ref. The detector is 

positioned in the far-field zone at f = f{0, 1,0}. All the lines consider single site addressing which is compared to the Mollow 
line (dotted black curve). In (c) and (d), small black dots indicate lattice sites, colored dots indicate the atoms considered for 
plots of the corresponding color in (a) respectively (b) and the faint orange circle indicates the FWHM of the probe laser. 



If the atoms in the configurations of Fig. HJc) and 
Fig. |4jd) are all driven by the same driving strength we 
find symmetric power spectra in all cases. In the case of 
two atoms the degree of asymmetry is exactly zero. In 
the cases of three or five atoms we find a degree of asym- 
metry which is negligible but not exactly zero. For two 
equally driven atoms the master equation is fully sym- 
metric under the exchange of these two atoms. This is 
not longer true in the case of three or five equally driven 
atoms as for example J 12 ^ J 23 because of the geometry 
dependence of the dipole-dipole interaction. However, if 



we choose J?" = J and P 1 " = T for all /i, u G {1, N} 
artificially, and drive all atoms with the same Rabi fre- 
quency, we find a spectrum for 2,3,4 or 5 atoms with a 
degree of asymmetry that is exactly zero. In the next 
section we discuss this observation further. 
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V. SOME CONCLUSIONS AND REMARKS 
ABOUT THE SYMMETRY OF THE SPECTRA 

Our results lead us to a remarkable first conclusion, 
regarding our particular system: 

1. In the case where the master equation and there- 
fore also the density matrix, is fully symmetric un- 
der the exchange of each possible pair of two-level 
atoms [i and v (i.e. under atomic permutation), the 
total spectrum of emission of the system is symmet- 
ric around the laser frequency. 

Moreover, putting together these results with many 
other examples one could think of [23], we further en- 
vision a second conclusion regarding any open quantum 
system in general: 

2. Provided that a single quantum system, QS, (such 
as a few-level system or harmonic oscillator under 
coherent or incoherent continuous excitation) ex- 
hibits a symmetric steady state power spectrum, 
then a number N of such QSs coupled to each 
other, in a way that the density matrix is fully 
symmetric under all possible permutations of these 
QSs, will also exhibit a symmetric total steady 
state power spectrum. Consequently all possible 
auto- and cross-correlation functions between pairs 
of QSs will be real and therefore experimentally 
observable. 

The second conclusion is a generalization of the first 
one. These statements are difficult to proof starting from 
the properties of the Liouvillian of an ensemble of coupled 
QSs. Here we only explore some directions that such a 
proof may take. 

The power spectrum for N QSs with associated oper- 
ator cr M , consists of a sum of N 2 contributions Sp, v (u), 
each given by 



Sft V (u}) oc 5? 



dr 



(17) 



with (cr+(r)cr~) ss = lim t _ >00 (<T+(t + r)a~(t)) and ui = 
lu — u>l- Each of these terms can be decomposed into a 
symmetric and an asymmetric part simply by separating 
the corresponding correlator into its real and imaginary 
parts: 



S MV (S>) = S%,(Q) + Sff(Q) = 

i 

dr $1{(<t+(t)<t~) ss } cos(wr) 



(18) 



In the case where the master equation, and therefore also 
the density matrix, is fully symmetric under QS permu- 
tation, we have S^^Cj) = S v ^{Cj). Consequently, the to- 
tal spectrum can be computed in terms of two different 



correlators only, for example 

S(Q) = NSn(u) + N(N - l)S 12 (u) • 



(19) 



In this case the spectrum is symmetric if and only if ev- 
ery correlator (er+ (t)<t~) ss is real. Then, the two-time 
operator cr+(r)cr~ with fx ^ v becomes an observable in 
the steady state as we stated in our second conclusion. 

However, the fact that the correlators become real 
when the spectrum is symmetric does not provide new 
information about the system, it is only a mathematical 
implication. The general question when an open quan- 
tum system should have a symmetric spectrum around 
some relevant frequency (the laser in the case of coher- 
ent excitation) is not a trivial one. First, this depends 
crucially on the nature of the excitation that is being de- 
tected, that is, the operators appearing in the two-time 
correlator. In our case generalized to N atoms, we re- 
fer to the collective Dicke operator = ^uLi a p but 
it could be any deexcitation operator in the system that 
corresponds to some physical entity. Each peak that ap- 
pears in the spectrum is related to the probability am- 
plitude to transit between two eigenstates of the system 
by emitting one of these quasiparticles. This means that 
both the dynamics of the dressed states and their quasi- 
particle component play a role. In order to make this 
link clearer, let us decompose the incoherent part of a 
spectrum into a sum of d 2 peaks poj , 



S(Q) 



1 



:E 
P =i 



L p j p /2 - K p (u - uj p ) 



( 7p /2) 2 + (w - u p f 



(20) 



with Up, 7 P (peak position and linewidth), L p and K p 
(Lorentzian and dispersive weights) all real parame- 
ters and d the dimension of the Hilbert space. Then, 
— (iujp+jp/2) are the eigenvalues of the Liouvillian in ma- 
trix form, L, which are either real (giving rise to a single 
Lorentzian peak at the center) or pairs of complex conju- 
gates (giving rise to a pair of sister peaks symmetrically 
placed around the center, with equal broadening) pH ]. 

A given pair of sister peaks (with iu> a + 7 Q /2 = 
— iup + 7/?/2) is symmetric if the complex weights are 
also conjugates, 



iK, 



P 



(21) 



These are computed from the eigenvectors of L includ- 
ing the steady state density matrix. They correspond 
exactly to the transition probability amplitude as men- 
tioned above. The balance between two amplitudes give 
rise to a pair of symmetric twin peaks. If this is the 
case for all sister peaks, the total spectrum is, of course, 
symmetric. 

In order to grasp all the physical sense of this balance 
condition, one would need to identify (or rather recon- 
struct) the eigenstates of the system under study. This 
is not an easy task, especially in presence of both dis- 
sipation and excitation [22l. [23j. In some limiting cases, 
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however, such as one strongly driven two-level system, 
it is possible 0, HE|- The eigenvectors of the full Li- 
ouvillian correspond in good approximation to the so- 
called dressed states, |±), obtained from diagonalizing 
the Hamiltonian part only. As in this regime the three 
peaks that form the Mollow structure are well sepa- 
rated, the interference part of the spectrum is negligi- 
ble (K p w 0) and the condition for symmetry of the 
two side peaks, L a = Lp, is completely equivalent to 
the so-called detailed balance between the dressed states: 

/0 ++ P( + _>_) = p--V( where p±± are the dressed 

state steady state populations and 7 , (±_ > z F ) the transition 
rates between them. On the other hand, if the excita- 
tion is weak and there is an overlap between the peaks 
of the spectrum the Hamiltonian, dressed states are not 
anymore the eigenstates of the full Liouvillian and their 
detailed balance is no longer a necessary condition for 
symmetry (in fact it breaks down out of resonance where 
the spectrum is still symmetric). 

A symmetric spectrum thus implies that the probabil- 
ity amplitudes of transitions between eigenstates are bal- 
anced. However, the challenging task of reconstructing 
such eigenstates makes it difficult to foresee and demon- 
strate when a system will exhibit a symmetric spectrum. 
Specially in a configuration where coupling strength, de- 
cay and excitation rates are of the same order of mag- 
nitude, as in this study. Putting Eq. (f2"T) in terms of 
the eigenvectors of the Liouvillian gives the mathemati- 
cal condition that they must fulfill so that the spectrum 
is symmetric. But this does not bring any further insight 
into the matter if one cannot identify which properties 
of the system cause the weights K a and L a of the corre- 
sponding Liouvillian L to fulfill this equation. 

On the other hand, our conclusions 1 and 2 are some- 
how intuitive and expected if one reasons on physical 
grounds, so let us end this Section with a plausible ex- 
planation for the symmetry of the spectrum and its break 
down in our particular case. 

In our configuration, the QS is simply a two-level atom 
driven by a laser, whose Mollow spectrum is indeed al- 
ways symmetric (in the absence of incoherent pumping, 
pure dephasing or other decoherence effects). When as- 
sembling N of such identical and identically driven QSs, 
new collective states are expected to appear. The driving 
in this case is restricted to the set of states |0) f-> erg |0) <-> 
... <-> (<7<t ) N \0) that form an (N + l)-level system. The 
remaining non-symmetric states are not driven at all but 
provide an effective decay channel for the (N + l)-level 
system. It is known from dressed state arguments that 
the total spectrum of a coherently driven (N + l)-level 
system is symmetric. Additional effective decay through 
the remaining non-symmetric states does not break the 
symmetry of the total spectrum, as we calculated in pre- 
vious sections. Now, if the atoms are not driven equally 
the laser does not solely drive the (N + l)-level system 
but also the transitions between non-symmetric states. 
This disrupts the dynamics of the (N + l)-level system 
and induces decoherence in the form of pure dephasing 



and an effective incoherent pump. Both elements are well 
known to break the symmetry in the spectra of coherently 
driven systems [H, [2g . Letting other atomic parameters 
be different, such as decay rates or detunings, has similar 
decoherent impact on the dynamics and the symmetry of 
the spectrum. We have checked that this is the case for 
systems consisting of up to five two- level atoms. 



VI. EXPERIMENTAL APPLICABILITY 

Our numerical results in section IIVI clearly show the 
signatures of single site addressing in resonance fluores- 
cence spectra. In summary it can be stated that the 
spectrum is symmetric around the laser frequency if the 
atoms are illuminated by a laser with equal strength. 
If, in contrast, only one atom is addressed the spectrum 
becomes asymmetric. Thus it is possible to measure a 
fluorescence spectrum and deduce information about the 
quality of an addressing scheme. 

In optical lattices a large number of ultra cold atoms 
are trapped. As soon as these atoms are illuminated 
by light the atoms heat up unavoidably. The detection 
methods presented in the experiments [8l-[lG| are all de- 
structive measurements in the sense that the atomic sam- 
ple is too hot after detecting the atoms. The single site 
addressing scheme presented in Ref. is designed in a 
way that the test for the failure or success of the address- 
ing scheme requires a destructive measurement. Hence 
the test and the usage of the addressing scheme have to 
be carried out in different atomic samples. 

Yet, if one makes a resonance fluorescence measure- 
ment as pictured in this work only a small fraction of 
atoms will be heated. Atoms positioned at the edges of 
a 2D lattice, for example, are well suited to test the ad- 
dressing laser. After the test of the addressing scheme 
the laser can be moved over the lattice to the desired po- 
sition. The test and the usage of the addressing scheme 
could be carried out in only one atomic sample. 



VII. CONCLUSION 

In this work we calculated and investigated the reso- 
nance fluorescence spectra of two-level atoms under the 
influence of local addressability. We implemented lo- 
cal addressing by means of a laser beam focused on a 
single atom in an atomic ensemble. A master equa- 
tion in a Markovian regime modeled the interaction be- 
tween the atoms and the surrounding electromagnetic 
vacuum. Due to small inter-atomic separations we had 
to account for a dipole-dipole interaction between atoms 
that is induced by the mutual exchange of photons and 
collective damping processes. With numerical calcula- 
tions we demonstrated that the output power spectrum 
of an atomic ensemble is asymmetric in the case of single 
atom addressing. We showed that this effect is gener- 
ated by the presence of the dipole-dipole interaction due 



9 



to photon exchange. Our results suggest that resonance 
fluorescence measurements could provide sensitive tests 
for the addressing of individual atoms in ID or 2D op- 
tical lattices. They thus also allow to predict emission 
spectra of ID optical lattices on the surface of optical 
nanofibres [29[. Our calculations are valid for any set 
of two- level quantum systems. The applicability of reso- 
nance fluorescence measurements as a test for single site 
addressability is therefore not restricted to neutral atoms 
in optical lattices but also applies to quantum dots [27j 
and colour centers in diamonds |28j ]. Furthermore, we 
provided some physical and intuitive explanation for the 
symmetry of the spectrum and its breakdown under sin- 



gle site addressing. We finally generalized our findings as 
follows: Provided that a single quantum system exhibits 
a symmetric steady state power spectrum, this property 
also holds for TV identical and identically coupled quan- 
tum systems. 
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